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An effective Coulomb gauge Hamiltonian, H e s, is used to calculate the light (uug), strange (ssg) 
and charmed (ccg) hybrid meson spectra. For the same two parameter H e g providing glueball 
masses consistent with lattice results and a good description of the observed u, d, s and c quark 
mesons, a large-scale variational treatment predicts the lightest hybrid has J PC = ++ and mass 
2100 MeV. The lightest exotic 1 _+ state is just above 2200 MeV, near the upper limit of lattice 
and Flux Tube predictions. These theoretical formulations all indicate the observed 1 _+ 7ri(1600) 
and, more clearly, 7ri(1400) are not hybrid states. The Coulomb gauge approach further predicts 
that in the strange and charmed sectors, respectively, the ground state hybrids have l" 1 with 
masses 2125 and 3830 MeV, while the first exotic 1" + states are at 2395 and 4020 MeV. Finally, 
using our hybrid wavefunctions, dimensional counting rules and the Franck-Condon principle, novel 
experimental signatures are presented to assist light and heavy hybrid meson searches. 



I. INTRODUCTION 

Following the "Eightfold Way" (Gell-Mann and 
Ne'eman), the "Quark Model" (Gell-Mann and Zweig), 
along with subsequent extensions, has generally ex- 
plained the observed hadronic spectrum. This is espe- 
cially true for heavy flavored mesons where it is now 
clear that higher order QCD corrections can be ignored or 
treated perturbatively. Even in the light sector, the phe- 
nomenological quark model works reasonably well. How- 
ever, the existence of hadrons with exotic quantum num- 
bers (i.e. J PC states not possible in qq or qqq systems) 
clearly reveals this model is not complete. Related, it is 
expected that there are exotic hadrons with conventional 
quantum numbers that also can not be described by the 
quark model, e.g., glueballs gg, hybrid mesons qqg and 
tetraquarks qqqq. 

Possible experimental evidence for a 1 h exotic state 
was first reported in 1988 [lj, but the situation was not 
clarified until several years later. Now it is believed that 
there exist two states with these quantum numbers below 
2 GeV: ^_l(1400) HOI and tti(1600) H0(note, a recent 
analysis 6] finds no evidence for either candidate) . There 
are also other reported hybrid candidates with J PC = 
0-+ Hi, 1— and 2"+ OSIIJIll. 

Theoretically, the structure of the tt\ states remains 
unclear. They could be hybrid or tetraquark mesons with 
most theoretical studies fl3lfl3 in ves tigating the former. 
Lattice gauge simulations |l5l ll6L UtI lisl Il9| predict the 
lightest hybrid meson is between 1.7 and 2.1 GeV and 
results from the Flux Tube model [2(], l2ll |22^ | also span 
much of this range. Only vintage Bag model |2^| calcu- 
lations yield a lower mass, between 1.3 and 1.8 GeV, but 
Ref. ;24J argues that the tt\ (1400) is not a hybrid. Table 
[IJlist predictions for the u/d, s and c 1 ^ hybrid mesons. 

In this work, we study qqg hybrid states using a field 
theoretical, relativistic many-body approach based upon 
an effective QCD Hamiltonian, _ff c ff, formulated in the 



Model [Reference] 


u/d hybrid 


s hybrid 


c hybrid 


Lattice QCD [15-19, 25-27] 
Flux Tube [20.,^!, 22} 
Bag Model [23j 


1.7 - 2.1 

1.8 - 2.1 
1.3 - 1.8 


1.9 
2.1 - 2.3 


4.2 - 4.4 
4.1 - 4.5 
3.9 



TABLE I: Published predicted 1 + masses, in GeV, for light, 
strange and charmed hybrid mesons. 



Coulomb gauge. This model successfully describes the 
meson spectrum [2|| and is also consistent [^3 with 
lattice glueball (and oddball) predictions. Using stan- 
dard bare current quark masses, it properly incorporates 
chiral symmetry, yet dynamically generates a constituent 
mass and spontaneous chiral symmetry breaking |3l| . 
Further, it provides a good description of the vacuum 
properties (quark and gluon condensates), and respects 
the global, internal symmetries of QCD, as well as the 
spatial Euclidean group, all within a minimal two pa- 
rameter theory. Our work also extends an earlier hybrid 
calculation 32] by including previously omitted terms in 
the Hamiltonian and by comprehensively predicting the 
light, strange and charmed hybrid meson spectra. 

This paper is organized into 8 sections. In Sections II 
and III the effective Hamiltonian is presented along with 
an improved hyperfine interaction which provides real- 
istic spin splittings in both light and heavy mesons |2^] 
and, for the first time, the rigorous non-abelian contribu- 
tions from the color magnetic fields. Section IV details 
the corresponding improved quark and gluon gap equa- 
tions and a variational formulation for the hybrid meson 
problem is developed in Sec. V. Calculations and new 
results are discussed in Sec. VI, while in Sec. VII we 
develop novel experimental signatures for observing hy- 
brid mesons having both conventional and exotic quan- 
tum numbers. Finally, we summarize results and conclu- 
sions in Sec. VIII. 
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II. EFFECTIVE HAMILTONIAN 



Confinement is described by a Cornell type potential, 



Our effective, quark-gluon Hamiltonian is an approx- 
imation to the exact Coulomb gauge QCD Hamiltonian 
|33|| and is given by (summation over repeated indices is 
used throughout this paper) 

H e g = H q + H g + H qg + He (1) 

H q = J dx& (x)[-ia ■ V + 0m]V(x) (2) 

H g = - J dx[n a (x) ■n a (x)+B a (x) -B a (x)} (3) 

H qg = gj dxJ a (x)-A a (x) (4) 

He = -\J dxdyp a (x)V(\x-y\)p a (y) . (5) 

Here g is the QCD coupling, is the quark field with 
current quark mass to, A a are the gluon fields satisfying 
the transverse gauge condition, V • A a = 0, a = 1, 2, ...8, 
11° are the conjugate fields and B° are the non-abelian 
magnetic fields 



B a = V x A a + ^gf abc A b x A c 



(6) 



The color densities, p a (x), and quark color currents, J a , 
are related to the fields by 

p a (x) = ¥{x)T a y(x) + f abc A b (x) -U c (x) (7) 
J a = ¥(x)aT a ^(x) , (8) 

where T a = 4j- and f abc are the SU3 color matrices and 
structure constants, respectively. 

The bare parton fields have the following normal mode 
expansions (bare quark spinors u,v, helicity, A = ±1, and 
color vectors ee=i,2,3) 



f dk 

* W " /<£? 

*e(*) - u x (k)b xc (k) + v x (-k)d[ c (-k) (10) 
dk 1 



A a (x) 



^c(k)e ikx e c (9) 
f UA (-fe)4 c (-fe) (10) 
[a a (fc) + SL ajt (-k)]e lk - x (11) 



(2tt)3 V2fc L 

n a (x) =-i |^y|[a a (fc)-a a t(_fe)]e lfe - (12) 

with the Coulomb gauge transverse condition, k a a (k) = 
(-l)"A M a%(fc) = 0. Here b xc (k), d xc (-k) and o°(fc) 
(p = 0, ±1) are the bare quark, anti-quark and gluon 
Fock operators, the latter satisfying the transverse com- 
mutation relations, 

K(fe),ajJ(fe')] = (2n) 3 5 ab 5 3 (k-k')D^(k) , (13) 



with 



D, p ,(k)=s^-(-ir^^ 



(14) 



V(r = \x-y\) = V c (r) 

V c (r) = -- 
r 

V L {r) = ar, 



V L (r) 



(15) 
(16) 

(17) 



where the string tension, a = 0.135 GeV 2 , and a s = 

2 

f- = 0.4 have been previously determined. The Fourier 
transform of V is denoted by V and in momentum 
space these potentials take the form Vl(|p|) = — 87rcr/p 4 , 
Vb(|p|) = — 47ra s /p 2 . For comparison and also to pro- 
vide hadronic structure sensitivity to potential form, we 
report predictions using a confining potential [34| having 
a renormalization improved short-ranged behavior. This 
potential was utilized in a previous meson study |2jfl and 
has the momentum space representation 



8.04 
7& 



l n -° l32 (Jl Tr+ 0.82) 



V(\ P \) = { ln°- a (^ + 1.41) P>m 9 (18) 



12.25mi 



P < m g 



The parameter m g sets the string tension and is related 
to dbym s ~ T/87TCT/12.25 w 600 MeV. 



III. HAMILTONIAN ff 2 CORRECTIONS 

As mentioned above, a previous hybrid application 
used this Hamiltonian but set the QCD coupling, g, to 
zero. This truncation eliminated the quark-gluon inter- 
action, J a ■ A", or "hyperfine" term, Eq. (4), and also 
the non- linear (non-abelian) component of the color mag- 
netic fields, Eqs. (3, 6). Now, both are included so that 
the non-confining part of the Hamiltonian is consistent 
to order g 2 . 



A. Hyperfine correction 

Following [§i|], the Hamiltonian term H qg containing 
the J a • A° operators is included using perturbation the- 
ory to second order in g. Then, integrating over the 
gluonic degrees of freedom yields an effective quark hy- 
perfine interaction with a J a • J a form. This contribution 
is represented by the Feynman diagrams in Fig. ^ 

The resulting transverse hyperfine interaction is 



1 



dxdyJ9(x)U ij (x,y)Jf(y), (19) 



V T = 



where the kernel reflects the transverse gauge 



Uij(x,y) = [ % - 



U(\x-y\). (20) 
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FIG. 1: The hyperfine correction entails the exchange of a 
gluon between q and q and qq annihilation. 




FIG. 2: The non-abelian correction with triple gluon vertices. 



For U we choose a modified Yukawa potential which in- 
corporates a dynamical mass, m„ = 600 MeV, for the 
exchanged gluon as explained in 29] . Fourier transform- 
ing to momentum space, this continuous potential takes 
the form 



8.04 



l n -0-62(JL 2 . +0 .82) 



U(p) 



P ln 08 (^- + 1.41) 
24.50 



P > m g 
P < m g 



(21) 



Non-abelian correction 



Similarly, the non-abelian components of the color 
magnetic fields, gf abc A b x A c , in the kinetic energy are 
also included perturbatively through (J a • A a )(B a • B a ). 
The resulting non-abelian interaction is represented by 
the Feynman diagram in Fig. Inland is given by 



dxdyr bc 6 t]k e a M{y){vf ) U kh {x,y))A b { X )A c m {x) + 



Vna = 4 _ 

J^(y)(V^ ) ^(x))?7 tt (x,y)^(x) + ^(y)(vf ) ^(x))4( :C )^(x,y)]. 



(22) 



r 



where C/y is the same kernel appearing in the hyperfine 
potential. 



IV. GAP EQUATION 

Having defined the model Hamiltonian, the next step is 
to calculate the ground state. Since we are free to expand 
the field operators in any complete basis, we follow the 
Bardeen-Cooper-Schriffer (BCS) method and perform a 
Bogoliubov-Valatin rotation, 



B\c(k) = cosy&Ac(fc) 



A sin — g? ac (— fc) 



Dxc(-k) 



cos -^-d\c(—k) + A sin -^-h 



2 " ACV 

AC 



(*0 



a a (k) = coshe fe a a (fe) + sinh 6 fe a aT (-fc) , (23) 

which transforms the bare particle operators a a , b\c and 
d\c into the dressed, quasi-particle operators a a , B\c 
and D\c, respectively. Now the fields are 

* c (fe) = U x (k)B xc (k) + V x (-k)D{ c (-k) 
dk 1 



A a (x) 

n a (x) 



(2tt)3 V2^ 

uJ k . 



[a a (k) + a a1f (-k)]e 



ih-x 



dk 



ih-x 



^-[ a -(k)-cc a \-k)}e^ 



where uj k = ke~ 2Gk . Note that the dressed quark ex- 
pansion remains functionally invariant with respect to 
the bare case since the quasi-particle spinors have the 
inverse rotation 



ZA(fc) 



Vx(-k) 



cos -^u\{k) - \ sin-^v\{-k) 



1 

71 



VI - sin^fc <j ■ k xx 
-v x (-k) + AsinyU A (fc) 



1 

71 



-VI =r sm0fc tr • k xx 
VI + sin (f> k xx 



(24) 



(25) 



Here the quark gap angle, 4> k = <j>(k), is related to the 
BCS angle 9 k by tan(4>k ~9k) = m/k. The quasi-particle 
(BCS) vacuum, defined by B xc \n) = D\c\£l) = a°\to) = 
0, is connected to the bare parton one, &ac|0) = ^acIO) = 

o o. by 

|^ fe ) = ^ / ^ Atan%bIc(fc)4c( - fe) |o> 



in 



gluon) 



f dk 

e J (^o 3 



■tanhe fc D^(fe)a»t(fe)a"I(-fe) 



|0) . 



The BCS vacuum, = \flq Uark ) <8> \Qgiuon), now con- 
tains quark and gluon condensates (correlated qq and 
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gg Cooper pairs). Performing a variational minimiza- 
tion of the vacuum expectation value of the Hamiltonian, 
5(n\H c ff\Q) = 0, independently with respect to <j) k and 
ujk, yields the mass gap equations for each sector 



HYBRID MESONS 



ks k - mcfe = 



dq 



[(SkCqX - S q Ck)V(\k — q\) 



3 J (2tt) 3 

2c k s q U{\k - q\) + 2c q s k W(\k - g|)](26) 



1 l A 

(jju = k 



3 f dq 



4 / (2tt) 3 



V(\k-q\)[l + x 2 



2l 1~ fc 



3 2 f dq 1 — x2 



A 9 J (2tt) 3 a,, 



(27) 



where 



^(1* - ,i) s u(\k - g|) ^!±^|(l±^) i (28) 

with Sfc = sin<f> k , c k = cos(f>k and x = kq. The last term 
in Eq. H27[l originates from the non-abelian component 
of the gluon kinetic energy. Dimensional analysis of the 
above integrals reveals that the first equation is UV finite 
for the linear potential since Vl(|p|) = — 87rer/p 4 , but 
not for the Coulomb potential Vc(|p|) = — 47ra s /p 2 . In 
Eq. H27fl there are both logarithmical and quadratical 
divergences in the UV region and an integration cutoff, 
A = 4 GeV, determined in previous studies is used. 

Once the current quark masses are fixed, the gap equa- 
tions can be solved numerically for the quark and gluon 
gap angles. Using \q) = B AC (fe)t|fi) and \g) = a a ^\Q,), 
the quark and gluon self-energies are respectively 



e k = (q\H cS \q) = ms k + kc k - 
^Ij0\3 [[skS q + c q c k x]V(\k - q\) 
+2s k s g U(\k - q\) + 2c q c k W(\k - q\)} (29) 
and, for fixed color index a (no sum), 



tlO\ - 



Sun' 



UJ q LU k 



2u k 
-D^(q) 



2 fJ 



dq 



(2tt) 3 2oj k u q 
[2D^{k)-D liV {k)D vv ,{q)D v , lll {k)\ , (30) 



both of which are infrared divergent in the presence of 
an infrared enhanced kernel. This is a welcomed feature 
of this approach, as colored states are removed from the 
spectrum. The infrared divergence cancels however in 
bound state equations for color singlet states leading to 
a physical spectrum of mesons and baryons. 



In previous publications |28|, |29j, |3fJ, |31| we have used 
this model to study the two-body meson and glueball 
systems by diagonalizing iJ c g using the Tamm-Dancoff 
and Random Phase approximations. We also made pre- 
dictions for three-body glueballs (oddballs) [3(j and pub- 
lished |32( a brief study of the three-body hybrid meson 
using a variational treatment. We now extend the latter 
and also provide more complete details of the variational 
calculation. 



A. Wavefunction ansatz and quantum numbers 

Following our initial study [3^ , we work in the hybrid 
center of momentum system and denote the momenta of 
the dressed quark, anti-quark and gluon by q, q and g, 
respectively. We then define q + = ^y 2 , q = q — q and 
note that g = q — q= —2q + . 

The color structure of a qqg hybrid is determined by 
SU(3) algebra 

(3® 3)0 8 = (801)08= (808)® (8 01) 

= 270 100 100 80 80 10 8. (31) 

Note for an overall color singlet the quarks must be in 
an octet state like the gluon. As discussed below, this 
leads to a repulsive qq interaction, confirmed by lattice 
at short range, which raises the mass of the hybrid meson. 
The hybrid wavefunction will therefore involve the color 
structure T^ iC J3^D s c< _o a ^ and has the general form 



I* 



JPC\ 



dg + dq_ JPC 
{2irf (2tt) 3 



^ lC X Cl (<z)< C2 (9K + (9)^} . (32) 



which is summed over color and angular momentum mag- 
netic sub-states. 

There are five angular momenta in this system, two or- 
bital, 1± (associated with q±) having z projections m±, 
and the 3 spins, S q = S q = 1/2 and S g = 1 with projec- 
tions Ai, A2 and fi, respectively. To form states with total 
angular momentum J, projection mj, we use the cou- 
pling scheme, S = S q + S q ,j = S g +l + ,L=j + l_,J = 
L + S. Then with the appropriate Clebsch-Gordan co- 
efficients, the hybrid wavefunction can be expressed in 
terms of a radial wavefunction F JPC '(<?+, <?_) and spher- 
ical harmonics, Y™ ± (q ± ), 

= F JPC {q + , q-)Yr + + (q + )Yr_- («_) X 

(-l)*- A2 (^,A 1 (-A 2 )|S'm s )(-l)^lZ +) (-/i)m + |im J -)x 
(jl-,mjm-\LrriL)(LS,rriLrns\JMj) . 

Since the intrinsic parity for a qq pair and a gluon are 
both —1, and the two orbital parities are (— and 
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1+ 


l- 
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Forbidden 
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1 






1" 




1 
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2~~ 




1 





1 








1 
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1-+ 


Forbidden 


1 
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1 


1 
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0" + 




1 
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1 
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1-+ 


Exotic 


1 





1 


1 


1 


2 
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2" + 




1 
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2 


2 


1 




+ 


1-+ 


Exotic 


1 





1 


2 


2 


2 




+ 


2"+ 




1 





1 


2 


2 


3 




+ 


3~+ 


Exotic 



TABLE II: Hybrid meson quantum numbers up to J = 3. 
Note exotic states and states forbidden by transversality: 
Z+ = 1 cannot couple to j = 0. 



(-1)'+, the total hybrid meson parity is 

P = (_i)(-i)(-i)' + (-i)'- = (-1)'++'-. (33) 

Finally, exchanging all additive quantum numbers, as 
required by charge conjugation, yields a (— l) l ~ +s factor 



from the space and spinor qq components which needs to 
be combined with the phase of the qqg composite color 
component. Although the gluon octet is not an eigen- 
state of C-parity, each gluon has a qq octet partner with 
opposite C-parity, resulting in a —1 contribution for the 
combined [[3 <£> 3]g ® 8]i system. Therefore the hybrid 
C-parity is 

C= (-l)(-l)'- +s - (-l) 1+ '- +s . (34) 

The extra —1 phase, as compared to a conventional qq 
meson having C-parity (— l) l ~ +s , is responsible for gen- 
crating exotic quantum numbers for certain hybrid states 
(e.g. J PC = 1 '"). Table Hi) lists quantum numbers for 
the model hybrid states for J up to 3. Note the ex- 
otic quantum number states and also states forbidden by 
the Coulomb gauge transversality condition (gluon or- 
bital 1+ = 1 can not couple with its spin to produce 
3=0). 

B. Variational equations of motion 

We now compute the hybrid mass, Mjpc , for each J PC 
with special interest focusing upon the exotic states. In 
terms of the above variational wavefunction, and upper 
bound for the mass is given by 

{^ PC \H ca \^) 

= M sel f + M m + M qg + M qgq + M ggg . (35) 

Here, the subscripts indicate the mass contribution from 
the self-energy of the three constituents, M se if, the qq in- 
teraction, Mqq, the qg and qg interactions, M qg , the 2 nd 
order correction from the qgq and qgq vertices, M qgq , and 
the 2 nd order correction from triple gluon vertices, M ggg . 
The three-body expectation value entails twelve dimen- 
sional integrals which can be reduced to nine dimensions 
by working in the center of momentum. The detailed 
expressions are 



M, 



self 



dqdq_jp C \ 



(2tt) 3 (2tt) 3 AlA2 ^ 



D vv > (g)(e q 



D llv {g)D ll . v .{9)e v q v 



(36) 



A I, 



\v{\q> - q\)Ul q U xW V{, q+ - q _ q ,Vx,- q + V(\q + g\H iq V^V{, q+ ._ q ,U xW 



(37) 
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M„ 



dq dq dq' 

;2tt) 3 (2tt) 3 (2tt) 3 



q+q+Ug'+q^JPCi 



^i^(q,mi^(Q',Q)D,'^(q + m^(Q + -Q)V(\Q ~ q'\)Ul q ,Ux x 



y/LOq+q l-O q > + q 

Uq+q+Uq+q > ^JPC\ {„ -\^JPC 
\JWq+q U q+q 



(38) 



M. 



qgq 



^mn(g' - g)Wj[ ig Q! m W Aiq 'V|, 9+ __ 9 ,Q! n VA 2? + f/ m „(g + q)Ul iq a m Vx 2 qV{^ q+q _ q ,a„U K q . 



(39) 



M 999 = 



rfq tfq dq' 



(2tt)3 (2tt)3 (2tt) 3 V o;_,,_ 5W , +5 
|vil/* h (g - g') A^C-g* - g) - D^(g)D k ^ (-q' - q 

+ i{q- q')iU lh (q' + q)D„ k (g)D kfi ,(-q' - q) 

- iU kh (q - q') \{q' + q)iD „i(g)D k ^ {—q' -q) + (q + q)iD, lk (g)D^ {-q' - q) 



(40) 



In the above expressions, e q , e q and are the quark, 
anti-quark and gluon self-energies, respectively, evalu- 
ated at the indicated momentum [g = —q — q). A picto- 
rial representation for each type of contribution is given 
by the Feynman diagrams in Fig. [3J 

The above expectation values are then computed 
variationally using the separable radial wavefunction, 
F(q+,q—) = /(<?+, ct+)f(c[-, ck_), having two variational 
parameters, a+ and a_. We investigated two functional 
forms for /; a gaussian and a scalable, numerical solu- 
tion from our two body meson studies. In general, the 




FIG. 3: Diagrams for (^ JPC \H cfi |* JPC ). 



gaussian radial wavefunction, 

/(g±,a±) = e- a ± > x ± = ^- (41) 

a± 

provided better results (lower variational mass) for s- 
wave states when compared to the numerical one. This 
was also true for p-wave orbital excitations, provided the 
gaussian was multiplied by x± corresponding to l± = 
1. All integrals were calculated using the Monte Carlo 
method with the adaptive sampling algorithm VEGAS 
|3fij . The integrals were evaluated several times with an 
increasing number of points until a weight-averaged re- 
sult converged. The hybrid mass error introduced by this 
procedure is about ±50 MeV. For each J PC hybrid state 
we optimized the variational parameters a + and a_ to 
produce the lowest mass. In terms of the string tension, 
their values fell in the ranges < a + < 1.7 y/a and 

l.OVcF < a_ < 3.3y/a. 

VI. RESULTS: HYBRID MESON SPECTRUM 

A. Light hybrid mesons 

For the light hybrid calculation we used m = 5 MeV 
|36j | for the u/d current quark mass. Results are listed in 
Table ITTT1 which shows the ground state is the ++ non- 
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(I)J PC 


Mjpc (MeV) 


Mjpc (MeV) 






no corrections 


with g 2 corrections 




(0)0++ 


2080 


2135 




(1)0++ 


2065 


2100 


Ground 


(o)i+- 


2135 


2140 




(1)1+- 


2135 


2140 




(0)2++ 


2340 


2335 




(1)2++ 


2180 


2170 




(0)1++ 


2415 


2470 




(1)1++ 


2110 


2170 




(o)i-+ 


2500 


2525 


Exotic 


(1)1-+ 


2205 


2220 


Exotic 


(o)o— 


2275 


2280 


Exotic 


(1)0— 


2280 


2285 


Exotic 


(0)1"+ 


2370 


2400 


Exotic * 


(1)1-+ 


2370 


2400 


Exotic * 


(0)1"+ 


2760 


2790 


Exotic 


(1)1-+ 


2570 


2600 


Exotic 


(0)3"+ 


3030 


3040 


Exotic 


(1)3~+ 


2910 


2915 


Exotic 



TABLE III: Spectrum of light hybrid meson states. Error 
~ ± 50 MeV. *Isospin degenerate states. 



(NA) and hyperfine corrections for several states. Gen- 
erally, both effects were small (except the hyperfine cor- 
rection for charmed quarks, see below), roughly of the 
same order as the overall 50 MeV Monte Carlo error. In 
particular, the NA correction entailed several terms with 
different signs which tended to cancel. 

Our model exotic spectrum (see Fig. 0J) spans almost 
a GeV, between 2.1 and about 3 GeV, and includes pre- 
dictions for J up to 3. There are no exotic J = 2 model 
states in this region since they require a d-wave or two 
p-waves, both involving much higher excitations. 

Finally, we comment on an interesting isospin splitting 
effect. From Fig. annihilation terms only contribute 
to the hybrid mass if the qq pair has quantum numbers 
consistent with the interaction. This is satisfied when 
Iqq = and S = 1. The annihilation diagrams can in- 
crease the 1 = hybrid states by several hundred MeV, 
as detailed in Table ITtTI In other cases, the l" 1 and one 
of the 1 h states should be isospin degenerate, as we 
compute to within the Monte Carlo error. On the other 
hand, the states 0++ and 0"" are not expected to be de- 
generate but, within the error, they are. It may be that 
the isospin splitting is not zero, but rather is smaller than 
the numerical error. Note the annihilation process for s- 
wave, isoscalar quarks in a triplet spin state is analogous 
to e~e + annihilation in the triplet state of positronium. 



exotic scalar, followed by the triplet 1+ , 2++ and 1++. 
The lightest hybrid mass is 2.1 GeV. 

For exotic states, as can be seen from Tabled at least 
one p-wave in q + or q is required. Because the qq inter- 
action is repulsive for quarks in a color octet state, the 
excitation energy is less for a 1+ (gluon orbital) excita- 
tion than a Z_ (qq orbital) excitation since the quarks are 
further separated and experience a larger repulsive linear 
force. 

The lightest exotic state is the / = 1,1 h , with mass 
2.22 GeV. This is slightly higher than the Flux Tube 
model and lattice QCD predicted masses for this state 
which were between 1.7 and 2.1 GeV (see Table 0. 

We studied the effects from including the non-abelian 



M(GeV 


jBxotica^ 


^exotica 


3.0- 


3~ + 

— r + 


3- + 


2.5 

2.0- 


- _r — r + 
— 2 ++ — r + 

— i- ~~ °~ 

_ — ++ 


— r 




uug (1=0) 


uug (1=1) 



FIG. 4: uug spectrum. 



B. Strange hybrid mesons 

Table llVl summarizes results obtained for the ssg (hid- 
den strangeness) hybrid mesons using a bare strange 
quark mass of 80 MeV. Now, the ground state is given 
by the non-exotic pseudovector state l" 1 , with a 2.125 
GeV mass, not at all reflecting the 75 MeV additional 
quark flavor mass contribution (the hybrid calculation is 
only sensitive to current quark masses above 200 MeV). 
Our prediction is in good agreement with the Flux Tube 
model and slightly above the only lattice prediction (see 
Table HJ). In the exotic sector, the lightest state is given 
by , with mass 2.3 GeV. Although there are also hy- 
brid states with explicit strangeness, e.g. sug, we do not 
show predictions since the effect from the s/u quark mass 
difference is small. 



C. Heavy hybrid mesons 

Table shows the results for the ccg (charmonium) 
hybrid mesons using a charmed quark mass of 1.0 GeV. 
The ground state is given, again, by the l" 1 state, with 
mass 3.83 GeV, while the lightest exotic hybrid lies at 
4.02 GeV. These numbers are in reasonable agreement 
with previous lattice and Flux Tube predictions, as listed 
in Table IJ 

Note that the correction introduced in the charmed 
case by the g 2 terms is roughly 500 to 600 MeV, signif- 
icantly higher than in the lighter hybrid systems where 
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J PC 


Mjpc (MeV) 
no corrections 


Mjpc (MeV) 
with g 2 corrections 




1+- 


2095 


2125 


Ground 


0++ 


2045 


2140 




2++ 


2290 


2315 




1++ 


2325 


2420 




1"+ 


2350 


2395 


Exotic 


0~~ 


2270 


2300 


Exotic 


1-+ 


2440 


2485 


Exotic 


1-+ 


2760 


2820 


Exotic 


3"+ 


2995 


3030 


Exotic 



TABLE IV: Spectrum of ssg states. Error ^ ± 50 MeV. 





Mjpc (MeV) 
no corrections 


Mjpc (MeV) 
with g 2 corrections 




1+- 


3310 


3830 


Ground 


0++ 


3295 


3945 




2++ 


3410 


3965 




1++ 


3450 


4100 




1"+ 


3545 


4020 


Exotic 


0"" 


3510 


4020 


Exotic 


1-+ 


3590 


4155 


Exotic 


1-+ 


3985 


4565 


Exotic 


3"+ 


4065 


4615 


Exotic 



D. Sensitivity to potential and parameters 

One of our key findings using the Cornell potential is 
that the mass of the lightest hybrid, especially the exotic 
1 '"j is above 2 GeV. Because of the ramifications of this 
result for exotic state searches, we have performed an 
interaction sensitivity study by varying both potential 
forms and parameters. 

We first varied the parameters in the Cornell poten- 
tial to obtain a lower bound for our predicted exotic hy- 
brid mass. Results are shown in Table I VII for different 
Coulomb potential parameters, a s , and string tensions, 
cr, found in the literature. For any combination of val- 
ues consistent with previous studies [2^, H^, |3(], 0OJ it 
was not possible to reduce the light hybrid mass to 1600 
MeV. In particular, we tried 0.0 < a s < 0.4 and 367 MeV 
< y/a < 424 MeV. Indeed, to obtain a hybrid mass as 
low as 1600 MeV required an unphysical y/a —262 MeV. 



potential/parameters 


I — 1 udg hybrid 


ccg hybrid 


Cornell [Eq. JT5J] 






s/a = 367 MeV, a s = 0.4 


2220 


4155 


Jd = 367 MeV, a s = 0.2 


2390 


4415 


^fd = 367 MeV, a s = 0.0 


2540 


4645 


^ = 424 MeV, a s =0.4 


2555 


4525 


Renormalized [Eq. I|18|l] 






m g = 526 MeV 


2705 


4730 


m g = 607 MeV 


3010 


5130 



TABLE V: Spectrum of ccg states. Error « ± 50 MeV. 



TABLE VI: Calculated light and charmonium hybrid 1 
masses, in GeV, using different interactions. 



the average corrections are 25 to 50 MeV. This large ef- 
fect arises from the hyperfine correction to the quark and 
anti-quark self-energies (see Eqs. (|29l I36[) ). which is en- 
hanced for heavier quark masses as discussed further in 
Ref. t 29j. Related, and as illustrated in Fig. the 
charmonium hybrid spectrum now has a slightly differ- 
ent level ordering from the lighter hybrid spectra. 



M(GeV) 
3.0 

2.5 

2.0 



exotica 



4.5 



4.0 



3.5 



exotica 



-0 ,1 



ssg ccg 

FIG. 5: Low lying ssg and ccg spectra. 



Table IVII also lists predictions for the confining po- 
tential given by Eq. I|18fl f° r values of the parameter 
m g = ^8ira/12.25 corresponding to the two different 
Cornell string tensions a but with the same current quark 
masses (m„ = 5 MeV, m c = 1 GeV). Note that this inter- 
action yields udg and ccg hybrids that are heavier than 
those given by the Cornell potential. Most significantly, 
this potential also predicts the lightest exotic hybrid has 
mass above 2 GeV. If we use m c = 0.85 GeV, which pro- 
vides a reasonable description of the charmonium spec- 
trum, the 1 ^ ccg mass decreases to 4815 MeV for m g — 
607 MeV. 



VII. SEARCHING FOR HYBRID MESONS 

Discovering exotic hadrons is a major goal motivating 
the Jefferson Lab 12 GeV upgrade and is also being ac- 
tively pursued by other collaborations and facilities, such 
as Babar, Belle, RHIC, etc. For low energy investigations 
of light quark exotic systems there is, unfortunately, no 
clean energy scale demarkation, since Aqcd governs the 
momentum distributions in light mesons and the strange 
quark mass is of the same order of magnitude. The ob- 
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vious detection strategy is therefore to perform statis- 
tically accurate cross sections measurements to extract 
partial wave amplitudes with explicitly exotic quantum 
numbers not accessible to ordinary qq states. However for 
(hidden) exotics with conventional meson quantum num- 
bers, it will be difficult to establish their nature. Note 
certain Flux Tube model pfl| and lattice predictions indi- 
cate that p-wave hybrid mesons prefer to decay to hadron 
pairs with one hadron also having a p-wave, rather than 
to two s-wave hadrons with a relative motion p-wave, e.g. 
rjhi in an s-wave as opposed to 7T7t in a p-wave. It will 
be interesting to check this prediction experimentally. 

More germane to the results of this paper are high 
energy experiments where novel tests can be conducted 
based on the scale separation provided by either the high 
beam energy or the heavy quark mass. This is discussed 
in the next two subsections. 



A. Application of dimensional counting rules 

Dimensional counting rules |3?| predict a power-law 
production cross section behavior for a given state at 
asymptotically high energies. They are based on the 
requirement that in forming a bound state with an 
energetic quark, the other partons must acquire very 
small relative momenta consistent with the production 
hadron's internal momentum distribution. This becomes 
highly unlikely in energetic collisions and therefore the 
production cross section falls as a power law, with the 
exponent increasing with increasing minimum number of 
constituents. For a hidden hybrid, the wavefunction has 
the Fock space expansion 



Px\cc) + /3 2 \ccg) + 



(42) 



where the quantum numbers are conventional but the 
first coefficient /?i is presumably small. Therefore the 
power-law behavior will reveal the second term in this 
series and permit distinguishing this exotic meson from 
ordinary cc charmonium (same for bottomonium). To be 
specific in the following discussion we focus on the re- 
cently discovered V>(4260) whose nature is currently un- 
der debate. The same remarks apply to any other hidden 
hybrid meson candidate. 

a. Inclusive production First consider the inclusive 
production reaction e~e + — > ip + X. The virtual photon 
fragments into a cc pair, each carrying half of the total 
center of mass momentum and therefore each has an en- 
ergy equal to the beam energy E beam (in a symmetric 
collider). The dimensional counting rules apply in the 
limit in which the produced particle's energy approaches 
its threshold value, that is, 



E. 



E, 



1 



beam 



(43) 



a minimum of one more gluon constituent in the leading 
wavefunction, (1 — x) 3 [38j. That is, as the tp is produced 
with more energy (this can be determined kinematically), 
the production cross section decreases linearly towards 
the kinematical endpoint where the ip has the maximum 
available energy. Note that the hybrid production cross 
section decreases even more rapidly as a cubic polyno- 
mial. With a sufficient number of events this can be 
documented experimentally. 

b. Exclusive production Bodwin, Braaten and Lee 
|39j have recently examined the reaction e~e + — > ipij}. 
Their study focused on the ground state J/ip, but their 
arguments also apply to the production of excited vector 
charmonia or a J/tp accompanied by a "0(4260). They 
studied rptp production as a function of the (small) vari- 
able ( m c is the charmed quark mass) 



E b cam 

and find in the limit r — ► 0, the differential cross section 
at fixed angle is constant. However, a straight-forward 
counting rule application predicts the production cross 
section is suppressed by a power of s or, equivalently, 
two powers of r, if one of the two produced hadrons is 
predominantly a hybrid meson, that is 



da 



dcos 9 
da 

dcos 9 



constant , for cc 
constant 

5 , for ccg 



(44) 
(45) 



In this limit, the power law behavior for a conventional 
charmonium state is 1 — x and for a hybrid state with 



Establishing this behavior only requires measuring dou- 
ble charmonium production at three sufficiently high en- 
ergies. For non-vector mesons both cross sections are 
further suppressed by additional powers of r, but the ex- 
tra r 2 signature always marks the presence of one more 
constituent, and therefore tags the hybrid meson. This 
argument also applies to the production of light hybrids 
in high-energy electron-positron colliders. 

To conclude this subsection, we note that the power- 
law predictions can be modified by QCD logarithms at 
very high energies and that in the strict limit 7-^0 
or x — + 0, any cc admixture in the wavefunction would 
dominate production. However for current, available low 
energies, the production will still be dominated by the 
hybrid component if (3\ << (32, i.e. if a predominantly 
pure hybrid state is found. 



B. Distinguishing the hybrid from charmonium 

In this subsection we propose a novel method applica- 
ble to the decay of heavy quarkonium that enables iden- 
tifying a new state as either a radial excitation of con- 
ventional charmonium or a hybrid state. The method is 
based upon two key points. 

First, we note that, due to the gluon mass gap scale, 
a conventional qq ground state (e.g. a well established 
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\|/(4s) 

Leading radial wavefunction squared 



C3 

O 

a 

>- 

s 



0.5 1 1.5 2 2.5 
p(GeV) 

FIG. 6: Probability density, \tp{p)\ 2 , for the 4s charmonium 
state. This is the relative cc quark momentum distribution. 

J/ip or T, etc.) is lighter than a qqg ground state hybrid 
with the same flavor. Indeed, the ground state hybrid 
mass is more comparable to a radially excited quarko- 
nium state. For example in our approach the i/ , (4s) and 
the ground state vector ccg state have similar masses. 
Now different eigenstates of a hermitian Hamiltonian 
are orthogonal with the nth radial excited state having 
n — 1 nodes. Therefore, even though the total energies 
(masses) are similar, the relative momentum distribution 
of the quarks in excited charmonium looks quite differ- 
ent from the quark momentum distribution in the ground 
state hybrid (see Figs. Eland[7J|. 

The second point involves the Franck-Condon (FC) 
principle widely used in molecular physics. Franck and 
Condon were the first to appreciate that molecular elec- 
tronic transitions proceed too rapidly for the much heav- 
ier nuclei to respond. The FC principle is applicable 
whenever there is a mass scale separation between differ- 
ent particles. In the context of quarkonium this means 
that the light fields (pions, gluons, etc.) quickly rear- 
range and the heavy quarks do not appreciably change 
their momentum distribution in the decay. Hence, the 
relative momentum between the decay products directly 
correlates with the quark momentum distribution in the 
parent quarkonium. Unfortunately in the simplest 2- 
body decays such as 

ip(ns) -> DD, T(ns) -> BB 

the FC constraint is not relevant since in the center of 
mass frame the momentum of the final products is fixed. 
This leads to smaller wavefunction overlaps suppressing 
the decay somewhat. However in 3-body decays such as 

ip(ns) -> DDtt, T(ns) -> BBir , 

the FC constraint applies. The first reaction can be em- 
ployed to study the recently discovered ?/;(4260). The 



Hybrid \|/(ls) 

cc relative wavefunction squared 




p(GeV) 

FIG. 7: Typical quark relative momentum distribution in a 
ccg hybrid state with a mass near 4300 MeV. The extra mass 
energy, relative to a cc J/ip, corresponds to gluon field exci- 
tations in collective models such as the Flux Tube approach 
or in the quasi-particle (gluon) mass gap in the constituent 
picture, but not in nodal radial excitation for the relative cc 
motion (compare to the radially excited charmonium distri- 
bution in Fig. |§J. 



second will be useful in an envisioned Belle collaboration 
measurement to establish whether this excited bottomo- 
nium state is the predicted quark model 5s meson state. 

Thus, we contend that the relative momentum distri- 
bution between the D and D mesons in the DDtt system 
mirrors the momentum distribution of the quarks in the 
parent ip meson. Since the hybrid ground state wave- 
function does not have a node, the resulting momentum 
distribution for the DD subsystem is also node-less and 
thus smoother than that for a conventional radially ex- 
cited charmonium. Multiplying by the relevant phase 
space distribution for this decay, yields the momentum 
distribution in Fig. |H1 that can be observed experimen- 
tally in the heavy-quark limit. The maximum is in the 
mid-momentum region where phase-space is larger and 
the wavefunction is near a local maximum. 

However, since quarks are not infinitely heavy, the FC 
signature is modified due to the recoil of the quarks in 
the D meson, yielding a different momentum distribu- 
tion. For example, taking a quark relative momentum 
between 150 and 200 MeV, one obtains the ccg momen- 
tum distribution illustrated in Fig. [5]for the ground state 
vector hybrid, and the smeared final state DD momen- 
tum distribution plotted in Fig. ^| for radially excited 
charmonium. As can be seen, even after smearing, there 
is still residual structure information adjacent to the cen- 
tral peak for radially excited charmonium that is remi- 
niscent of its parent charmonium wavefunction behavior, 
in sharp contrast to the smooth, bell-shaped hybrid dis- 
tribution. 
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Y(4s) 

Squared wf times phase space 




0.5 1 1.5 2 

p(GeV) 

FIG. 8: Momentum distribution of Fig. |S| multiplied by 
the phase space for the decay ^>(4260) — > DDtt. This is the 
probability density for finding a DDn state with relative DD 
momentum p according to the Franck-Condon principle in the 
heavy-quark limit. This signature will be more robust for the 
related bottomonium process T(5s) -— > BBn. 



Hybrid y(ls) 

DD relative momentum distribution in DD7t 




p(GeV) 

FIG. 9: The relative momentum distribution of the DD pair 
in the DDtt final state for a charmed hybrid meson with mo- 
mentum distribution given in Fig. [7] The distribution of the 
final products has a smooth bell shape in sharp contrast to 
the radially excited quarkonium distribution in Fig. I1UI 



We therefore advocate analyzing the DD and BB rel- 
ative momentum distributions in DDtt and BBn decays 
of highly excited quarkonia. Additional final state pions 
or other light particles do not alter our arguments (but 
restrict somewhat the available phase space) , so there are 
several other possible final state channels to search. 



Y(4s) 

Squared wf times phase space 




p(GeV) 

FIG. 10: The momentum distribution in Fig. [§] (solid line) 
and the distribution (dots) obtained averaging over a 150 MeV 
spread to approximate the actual quark momentum distribu- 
tion in the D meson. The experimental signal, firmer for 
bottomonium than charmonium, is a central peak with two 
shoulders for the relative momentum distribution of the BB 
(or DD) pair in the BBn (or DDtt) final state. These adja- 
cent enhancements identify the resonance as a radially excited 
quarkonium state as opposed to a hybrid meson. 



VIII. DISCUSSION AND CONCLUSIONS 

Our key model predictions are that the lightest hy- 
brid mass is 2.1 GeV, with the lightest 1 h exotic state 
slightly above 2.2 GeV. Lattice and Flux Tube calcula- 
tions yield a mass of at least 1.8 GeV for the 1 '", ex- 
cept Ref. 0] which predicts a mass around 1.7 GeV. 
Thus, with exception to this last work, the compos- 
ite model analyses appear to preclude the possibility of 
the reported 1 h exotica, 7Ti(1400) and 7Ti(1600), be- 
ing hybrid mesons. If this is correct, one should inves- 
tigate other structures for those two hadrons, such as 
tetraquark molecules, with both qq pairs in color singlets, 
or tetraquark atoms, where quark pairs are in intermedi- 
ate non-singlet color states and we are currently applying 
our model to these systems. However, if the lattice 1.7 
GeV 1 ^ prediction is robust, one can not yet exclude 
the observed 7Ti(1600) from being a hybrid, but this still 
does not explain the structure of the 7Ti(1400). It would 
be very useful to have other lattice measurements, using 
the same techniques as Ref. ^8|, to confirm or reject 
this result. Related, we have also varied our model pa- 
rameters and potential forms to obtain a lower bound for 
our predicted exotic hybrid mass which is clearly above 
2 GeV. 

Regarding isospin splitting, our results show an en- 
hanced splitting from g 2 corrections. For the ++ hybrid, 
the corrections increased the splitting from 15 to 35 MeV, 
the maximum increase in the light hybrid spectrum. 
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In the strange sector, we predict the lightest non-exotic 
hybrid mass is 2.125 GeV, while the lightest exotic mass 
is 2.30 GeV. These values compare reasonably well with 
Flux Tube [S3, 123, IS^I and, slightly lighter, lattice [Is| 
results. For the charmed sector, our predictions of 3.83 
GeV for the lightest hybrid and 4.02 GeV for the low- 
est exotic are also in good agreement with several other 
lattice and Flux Tube studies (see Table [Q). 

As mentioned above, the different g 2 corrections pro- 
duced an overall small effect, about the same order as 
the Monte Carlo error. However, the hyperfine correc- 
tion becomes important for heavier quark mass. In the 
charmed sector, this correction added about 500 to 600 
MeV to the hybrid mass. Lastly, note the level ordering 
of the exotic isoscalar uug and ssg spectra are the same, 
0~~, 1 1 ^ 1 h and 3 h , but slightly different than 
the exotic ccg system, where the and lowest 1 <~ 
are degenerate. This is a consequence of the enhanced 
charmonium self-energy from the hyperfine interaction. 

Finally, we discussed both low and high energy scenar- 
ios for observing hybrid mesons. For low energy studies 
involving light quark systems, dimensional counting rules 
predict specific power-law behaviors for distinguishing 
between production of conventional and hybrid mesons. 
For high energy investigations of heavy quark systems, 



the Frank-Condon principle provides a useful constraint 
on the final state momentum distributions which should 
assist experimentalist in identifying heavy hybrid sys- 
tems. 

In summary, lattice, Flux Tube and our H e g many- 
body approach all predict similar hybrid spectra and 
that the lightest 1 *~ exotic hybrid meson mass is near 
2 GeV. This composite model agreement indicates that 
the 7Ti(1600) is not a hybrid meson but has an alterna- 
tive structure. If true and if the 7Ti(1600) exists, it is 
more likely a tetraquark system, either a (qq) (qq) meson 
molecule or an exotic qqqq atom. Future work will apply 
our model to light and heavy tetraquark systems includ- 
ing mixing with hybrid and conventional meson states. 
Three-body forces ^3] will also be examined. 
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